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ABSTRACT 

Unsteady activity of 7-ray burst sources leads to internal shocks in their emergent 
relativistic wind. We study the emission spectra from such shocks, assuming that they 
produce a power-law distribution of relativistic electrons and posses strong magnetic 
fields. The synchrotron radiation emitted by the accelerated electrons is Compton 
up-scattered multiple times by the same electrons. A substantial component of the 
scattered photons acquires high energies and produces e + e~ pairs. The pairs transfer 
back their kinetic energy to the radiation through Compton scattering. The generic 
spectral signature from pair creation and multiple Compton scattering is highly 
sensitive to the radius at which the shock dissipation takes place and to the Lorentz 
factor of the wind. The entire emission spectrum extends over a wide range of photon 
energies, from the optical regime up to TeV energies. For reasonable values of the wind 
parameters, the calculated spectrum is found to be in good agreement with the burst 
spectra observed by BATSE. 
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1. Introduction 

The detection of Fell and Mgll absorption lines at a redshift of z = 0.835 in the optical 
spectrum of GB970508 (Metzger et al. 1997), provided the first confirmation that 7-ray bursts 
(GRBs) originate at cosmological distances. Most of the qualitative properties of cosmological 
GRBs are explained by the fireball model (see e.g., Goodman 1986; Paczyhski 1986; Meszaros <fe 
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Rees 1993 [MR]). In this model, a compact (~ 10 6 ~ 7 cm) source releases an energy of ~ 10 52 erg 
over a duration <^ 10 2 seconds with a negligible baryonic contamination (<; 1O~ 5 M ). Unsteady 
activity of the source results in a wind composed of many thin layers (fireball shells) of varying 
energy and baryonic mass. Within each shell the high energy-density at the source results in an 
optically thick e + e~-pair plasma that expands and accelerates to relativistic speeds. After an 
initial acceleration phase, the radiation and thermal energy of the fireball plasma is converted into 
the kinetic energy associated with the radial motion of the protons. Collisions between the shells 
can convert part of that kinetic energy into radiation and yield the primary GRB via synchrotron 
emission and inverse-Compton scattering (Paczyhski & Xu 1994; Rees & Meszaros 1994 [RM]; Sari 
& Piran 1997 [SP]). As the wind continues to expand, it impinges on the surrounding medium and 
eventually drives a relativistic blastwave in it, which heats fresh gas and accelerates electrons to 
relativistic speeds, thus producing the delayed afterglow radiation observed on time scales of hours 
to months (van Paradijs et al. 1997; Bond 1997; Djorgovski et al. 1997; Mignoli et al. 1997; Frail 
et al. 1997) via synchrotron emission (Wijers, Rees, & Meszaros 1997;Waxman 1997a,b; Vietri 
1997a,b) . 

The primary GRB emission is more likely caused by internal shocks than the external shocks 
(MR), since they can occur closer to the source and thus account for the rapid variability observed 
in many bursts (RM; SP). Unsteady activity of the central source naturally results in faster shells 
overtaking slower ones in front of them, and hence in energy dissipation by internal shocks. The 
complex temporal structure observed in GRBs then reflects the activity-history of their sources 
(SP; Kobayashi, Piran, & Sari 1997 [KPS]). 

It is often assumed that behind internal shocks, electrons are Fermi accelerated with a near 
equipartition energy density and magnetic fields acquire nearly equipartition strength. The 
electrons cool by synchrotron emission and inverse-Compton (IC) scattering off the synchrotron 
photons. Under typical conditions, the time scale for IC scattering is shorter than the synchrotron 
cooling time. Multiple scattering of the photons boosts a significant fraction of the radiation 
energy to frequencies above the e + e~-pair creation threshold. The pairs produced in this process 
are also relativistic and cool rapidly by IC scattering. Since the annihilation time of these pairs 
is longer than the hydrodynamic time in the comoving frame, they survive in the wind for a long 
time. Although the creation of pairs and their subsequent cooling is likely to leave noticeable 
imprints on the emergent radiation spectrum, it has not been analyzed before in the GRB 
literature. Since the photon and electron densities decline rapidly with radius, the strength of 
these signatures can serve as a probe of the radius at which the internal shocks occur. 

The GRB spectrum should also depend on the level of baryonic contamination in the wind. 
The extreme limit of pure energy release with no baryons was ruled out in the past based on the 
prediction that a point explosion of this type would lead to a roughly thermal spectrum (Goodman 
1986). One can place a lower limit on the baryonic mass in the fireball shells by requiring that 
internal shocks should occur before an external shock does (this limit depends on the ambient 
medium density). On the other hand, an upper limit can be placed based on the variability time 
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scale of the source and the condition that the shells be optically thin at the radius where internal 
shocks occur. 

In this Letter we study in detail the emergent spectra from the collision of two fireball shells^]. 
In particular, we quantify the significance of the radiation processes which were previously ignored 
in the literature, such as multiple Compton scattering, e + e~ creation, and their subsequent 
cooling. We use the collision kernels, reaction rates, and the computational techniques given in 
Pilla & Shaham (1997 [PS]). More details about this calculation and an elaborate study of the 
spectral characteristics of relativistic shocks will be included in a subsequent publication (Pilla 
&; Loeb 1997). In §2 we describe our model and specify the physical conditions in the emission 
region. In §3 we outline the relevant radiation processes and compute the model spectra. Finally, 
§4 summarizes the main implications of this work. 



2. Physical Properties of Internal Shocks 

The typical fireball dynamics (Piran, Shemi, & Narayan 1993, [PSN]; Meszaros, Laguna, & 
Rees 1993; Sari, Narayan, h Piran 1996 [SNP]; SP; and KPS) can be illustrated by considering 
a single shell of total energy E = lO 51 ^! erg, rest mass M = 10 27 M27 g, and initial radius 
ro = 10 7 r7 cm. After a brief acceleration phase, the Lorentz factor of the shell reaches a constant 
value r at an observer-frame radius r m ~ Tr$ (the protons are taken to be non-relativistic in 
the comoving frame, before the collision of shells). The energy of the shell is predominantly 
kinetic beyond this stage. Outside the radius r ~ T eX p — 2r 2 ro (all radii in the present analysis 
are measured in the observer's frame), the comoving width of the shell increases linearly with 
radius (PSN). The comoving proton density scales cis ti oc t for shell radii v v eX p and as 
n oc r~ 3 for r > r exp . The Thomson optical depth of the shell is r ~ (rt/r) 2 for r > r exp , where 
rt ~ (Mo"T/47rm p ) 1 ^ 2 ~ 6.3 x W^M^ 2 cm is the radius at which the shell becomes optically 
thin to Thomson scattering by its own electrons. Here <jt is Thomson cross section and m p is 
the proton rest mass. For r > r exp , the comoving width of the shell is 5 ~ r/2T. Two photons 
which are emitted with a proper-time difference ~ S/c, reach the observer with a time separation 
~ 5/2Tc ~ r/4r 2 c (SP). Assuming that the radiative cooling time is much shorter than the light 
transit time through the system, one finds that the observed width of the radiation pulse t p and 
the radius of the emission region r e are related through r e ~ 2T 2 ct p , which in turn can be used to 
constrain the values of 5, r, and n from observations. The emission spectrum will be nonthermal 
only if r e S> Tt- 

Now consider a wind of total duration td, composed of many thin fireball shells of thickness 



*In principle, each individual peak in the burst light-curve might correspond to a collision of a pair of shells 
(KPS) . The time-averaged flux of the entire burst is then an energy-weighted sum of the contributions from individual 
collisions. Here we provide the time-averaged spectra for the collision of two shells, whereas the generalization to the 
entire burst (with multiple peaks) is straightforward. 
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5 <C ctd (6 may change from one shell to another). We assume that for r > r exp there are regions 
in the wind where slower shells precede the faster ones, i.e. dT(r)/dr < 0, where T(r) is the 
local Lorentz factor. The spatial extent of the wind is typically ~ ctj ^> r$ (in fact, 5 S> r$ also 
holds in general). It was shown by Waxman & Piran (1994) that under these conditions the 
wind layers are susceptible to Rayleigh- Taylor instability because a rarefied fluid shell is pushing 
against a denser one. The resulting turbulent mixing will complicate the shock structure and 
deform it away from a simple planar geometry. Merging of a rarefied shell with a denser one might 
therefore be accompanied by the formation of "fingers" perpendicular to the shell walls, similar 
to the non-relativistic shock structure in supernova remnants (e.g., Jun & Norman 1996). The 
combined shells would then break into bubbles of different sizes, and the energy dissipation would 
take place near the bubble walls, due to collisions among them or instabilities on their surfaces 
(Kamionkowski & Freese 1992). 

By assuming that the shells as well as the shock fronts remain planar, KPS had found that 
the dissipation efficiency of internal shocks might obtain high values (~ 50%) for reasonable wind 
parameters. We assume that similar efficiencies are achieved in the case of unstable mixing. 
Since the combined area of the bubble walls greatly exceeds that of a planar shock, the electron 
acceleration efficiency in the present case is likely to be higher. However the temporal and spectral 
characteristics of the bursts might be different in the two cases. For a planar shock the accelerated 
electrons populate a thin layer around the shock front since their cooling time is much shorter 
than the transit time of sound waves across the shells. However, in the unstable mixing case 
we assume, to a first approximation, that the energy dissipation takes place nearly uniformly 
and simultaneously throughout the entire volume of the emission region. For the purpose of 
estimating the physical conditions involved, we take E » 10 51 erg, M » 10 27 g, r e rs 10 14 cm, and 
a bulk Lorentz factor of the emission region of V ~ 400. These values yield n ~ 3.1 x 10 10 cm~ 3 , 
r«3x 10~ 3 , a comoving width of the post-shock shell A « 1.4 x 10 11 cm, and an average Lorentz 
factor of the protons in the comoving frame (i.e., after shock heating) of j p s=s 3. 



3. Radiation Mechanisms and Spectra 

Acceleration of electrons to relativistic energies and the presence of strong magnetic fields 
are essential for converting the energy dissipated by the shock waves into radiation. Since the 
physics of neither of these processes is well understood, we parameterize the corresponding 
energy densities in units of their equipartition values. A magnetic equipartition parameter 
Qb = -B 2 /87rnm p c 2 (7 p — 1), corresponds to a field strength B 1.9[CBnio(7 p — 1)] 1 / 2 10 4 G, where 
n w = n/(10 10 cm -3 ) (all quantities are comoving, unless stated otherwise). We assume that the 
electrons are accelerated throughout the emission region. The amount of energy transferred from 
protons to electrons is uncertain, and we define the acceleration efficiency C, e in such a way that 
the average Lorentz factor of the electrons immediately after they are heated (via Fermi-type 
acceleration) is 7 = Ce(7 p — l)m p /m e , where m e is the electron rest mass. Because the Coulomb 
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collision time ~ 1/ cnax In A ~ A/er In A » A/c, we can safely ignore collisional relaxation in our 
analysis (here In A ~ 30 is the Coulomb logarithm) . 

Collisionless acceleration of electrons can be efficient if they are tightly coupled to the 
protons and the magnetic field by means of plasma waves (Kirk 1994). The typical Alfven speed 
in the plasma is va = B / (Aimm p ) 1 / 2 = min(l, 1.4^ 2 )c. The Larmor radius of an electron 
with a Lorentz factor 7 is rL, e ~ 0.l7(Csraio) -1 ^ 2 cm, and that of a proton of equal Lorentz 
factor, r^ iP , is larger by a factor m p /m e . The corresponding acceleration time scales (e.g., 
Hillas 1984) are, t acCje ~ cr LiC /v A ~ 1.7 x 10~ 13 7(C|raio)~ 1/2 sec and t acc , p ~ (m p /m e )t acc>e , 
respectively. Synchrotron losses limit the maximum value of the electron Lorentz factor to 
lmax,e ~ 3.2 x 10 5 (nf £B) _1//4 . The minimum value is determined by 7 and the shape of the 
electron distribution. We assume that the fraction of electrons per unit Lorentz factor 7 has the 
form F e (j) = (p - I)7m7n7~ p for <y min < 7 < ^ max . Thus "f min = (p - 2)7/(3? - 1), where 7 is the 
average Lorentz factor at any given time. The energy density in electrons immediately after their 
acceleration is uq = nm e c 2 ^fQ. In our model the radiation time scale is much shorter than the 
hydrodynamic expansion time scale in the comoving frame (which leads to a radiative efficiency of 
nearly 100%) and the radiation density at the end of electron cooling is therefore no- We compute 
the spectra for Qb = 0.1, Ce = 0.3, and p = 3.5. 

Electron Cooling and Radiation Spectrum 

We assume that the radiation energy density is initially small, and hence the electrons start 
losing their energy via synchrotron emission. As the energy density of the emitted radiation builds 
up, cooling via IC scattering becomes important as well. The typical time scale for synchrotron 
or IC losses is t c ~ q/ XcnaT'J 2 , where A = Cfi/Ce in the synchrotron case and A = u^/uq in the 
IC case. Here u-y is the radiation energy density and q is a dimensionless constant of order unity. 
For typical GRB conditions t acc ^ e <C t c <C to, where to ~ A/c is the light transit time through 
the system. Therefore electron cooling takes place only locally. We assume that the acceleration 
proceeds throughout the cooling phase in such a way that it maintains a steady power-law 
distribution with a constant index p, while 7 declines due to radiative losses. If Qb *C Ce, the 
synchrotron cooling time is long and multiple IC scatterings become important. Each scattering 
in the Thomson regime increases the photon energy by a factor 4(7 2 )/3 (Loeb, McKee, & Lahav 
1991), so that some of the photons are eventually boosted into the Klein-Nishina (KN) regime. As 
7 decreases, the up-scattered part of the spectrum spreads over many decades in frequency. The 
electron energy density u e = nm e c 2 ^f changes at a rate 

dt ~ \ dt ) sy + \ dt ) jc 1 U 

where {du e /dt) sy = —4ub/3T, ub = B 2 /8ir, 1/T = cnor(7 2 )o, and (7 2 )o is the average value 
of 7 2 over the initial electron distribution. The IC cooling rate, (du e /dt) IC , is derived below. 
The characteristic cooling time t c is a few times T under typical conditions. If u^it) is the 
instantaneous energy density of radiation, then energy conservation implies u 7 (i) + u e (t) = uq 
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(which is constant), and by assumption u 7 (0) = 0. The instantaneous radiation spectrum can be 
characterized by $(e,t) so that ii 7 (i)$(e, t)de = v(e,t)de is the fraction of the radiation energy 
density within an interval de around e; here e = hv/m e c 2 , v is the photon frequency, and h is 
Planck's constant. 



The spectral evolution rate (see Pilla & Loeb 1997) is derived from the equation 







du e 



uy(t)— *(e,t) - §(e,t) — e - 



dt 



dt 



9 , : 



*,</ 



ic 



(2) 



where we have used the fact that du^/dt = —du e /dt. The first term on the right hand side is the 
synchrotron emissivity from relativistic electrons (Rybicki & Lightman 1979) which depends on 
F e (7,t); the second term is 
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where fj, is the cosine of the scattering angle, 7; = max(7 m j n , e), £ = 7/(7 — e) = 1 + 27(1 — /i)ei, 
£1 = e/27(l — /x)(7 — e), and ukn = 3ot(C 2 — 2^/3 + l)/4£ 3 is the Klein-Nishina cross-section. By 
integrating both sides of equation (||) over all values of e, we obtain (duj/dt) IC on the left hand 
side, whereas on the right hand side we use the relation \de/dsi \ = 2(1 — /i)7 2 /C 2 to convert the 
integral to be over s\ and integrate over all values of 7, £1, and fi. Finally, we divide the result by 
2 to avoid double counting of each scattering event, and obtain 
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are the limiting energies of photons in the plasma. In the Thomson regime £ = 1 
cjt, and one finds 
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in agreement with the well-known result for 7 > 1 (Loeb et al. 1991). The coupled equations 
(|l])-(Q) are solved numerically for the radiation spectrum at the end of the cooling process (i.e., 
when u-y — » uq). Cooling ends after a relatively short time, <J 10T <C £q, for the parameters of 
interest here, and the photons decouple from the electrons subsequently. An example for the time 
evolution of the spectrum is shown in Figure la. It is evident that the radiation density above the 
e + e~-pair creation threshold is substantial. 

The spectral evolution rate due to pair creation is described by 
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where the integration is over the range — 1 ^ fx ^ 1 3>nd. £min < £ < £max subject to the condition 
ee'(l — /i) > 2. Here <r± is the pair creation cross section (PS) and 1/T" = cnax^Q- In all examples 
considered in this Letter, we find that t c 2i(e,t)/T' <C 1 for all relevant photon energies. Therefore, 
pair creation comes into play after the original electrons cool and decouple. The total energy loss 
per unit time and volume due to pair creation is obtained by integrating equation (|6|) over all 
values of e, 

The rate at which pairs are created per unit volume is 

/ dn\ n f £ ™<™ n ~ , , 

\*h = vL <^(M)^(t), < 8 > 

where Z^e, i) is same as Xi(e, t) except that e/70 is absent in its integrand. At the beginning 
of pair creation, the energy density of radiation is uo. The average Lorentz factor of the newly 
created pairs is 7±(i) ~ 7 y Ii(t)/l2(t), since u^/n pa 7 . Because the pairs are born relativistic, 
they transfer almost all their energy to the radiation via IC scattering (and also synchrotron 
emission if Qb ~ Ce) on a time scale of order a few T'. Similar cascades of pair creation and cooling 
occur in active galactic nuclei (Svensson 1987). Throughout the pair cascade and cooling process 
the radiation spectrum evolves continuously but its total energy density remains nearly constant 
at a value close to u$. Figure lb shows an example of this evolution up to a time ~ to- Because 
the hydrodynamic time scale in the comoving frame thyd ~ r/cT pa 2to, pair processes can not 
continue to operate for times much longer than to because of the decline in the densities due to 
the expansion of the fireball. 

Spectral Predictions of the Model 

We solved the coupled equations given above using the methods described by PS and 
obtained the model spectra for E pa 1.1 x 10 51 erg, 7 p = 3, Cb = 0.1, Ce = 0.3, and p = 3.5. 
We assume a source redshift z s pa 1 and include its effect on photon energy and flux. The total 
observable fluence for the above parameters is ~ 3.1 x 10~ 7 erg/cm 2 . 

Figure 1 shows the entire spectral evolution (in the comoving frame) due to the synchrotron 
and IC cooling and pair cascade processes. We take M w 10 27 g which yields T ps E/*f p Mc 2 w 400 
( which are consistent with empirical constraints, e.g. Woods & Loeb 1995); and a dissipation (or 
shell collision) radius of 3.3 x 10 13 cm. Note that the assumed values for E and M correspond 
only to a single emitting shell; they should obviously be higher for the entire burst. For the above 
parameters, we get A w 4.1 x 10 10 cm, n pa 1.1 x 10 12 enr 3 , r pa 3 x 10~ 2 , and B pa 6.3 x 10 4 G. 
The characteristic time scale for synchrotron and IC cooling is T = 1.6 x 10 _5 to- We find that 
after a time t = t c pa 16T, only 10% of the initial energy remains with the electrons. Figure la 
illustrates that as a result of electron cooling, the location of the synchrotron peak shifts to lower 
energies and also becomes broader. Later on, another peak develops at a much higher energy due 
to IC losses, with KN suppression in the very high-energy tail. Figure lb shows the evolution 
much later due to the pair-cascade process, which depletes the population of high-energy photons. 
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The IC loss by these pairs produces a power law tail at yet higher energies. The low-energy 
synchrotron part of the spectrum is almost unaltered throughout this phase. The computation 
was stopped at t « to at which time most of the photons leave the system. 

Figure 2 shows the dependence of the emergent spectrum on the dissipation radius for a fixed 
value of r (panel [a]) and dependence on T for a fixed value of the dissipation radius (panel [b]). 
Panel (a) shows that the location of synchrotron peak moves up in energy as the radius decreases 
(with an opposite trend for the IC peak); hence the overall extent of the emission spectrum gets 
narrower at smaller radii. Since B oc n 1 / 2 cx r~ 3 / 2 and r oc r~ 2 , as r decreases the synchrotron 
frequency increases and the pair-production depletion of high-energy photons is enhanced. One of 
the model spectra in Figure (2a) is compared with the BATSE data through the empirical formula 
for the time-integrated flux given by Band et al. (1993 [B93]). In that case our model predicts 
~ 50% more time-integrated flux at 10 keV than observed, which could be due to an oversimplified 
form of the electron distribution function which we have adopted. In our notation, the number of 
photons per unit interval of e is dN/de oc e~ 1 ^(e). Therefore, Band's formula (B93) reads 

e _1 $(£) = c\E~ a exp(— e/eo), if e < ((3 — a)eo 

= C2E~ lS , otherwise. (9) 

The constants are fixed by the requirement that this function be continuous at e = (/? — a)so 
and the normalization Jq 1 de$(e) = 1. Note that we have altered the signs of both indices a and 
(3 relative to the convention of B93. The observed values of the parameters (cf. Table 4 in B93) 
are in the range 0.3 < a < 1.5 (although in a few cases a is zero or negative), 1.6 < (3 < 5, and 
15 < eoirieC 2 (keV) < 3 x 10 3 . The majority of the bursts seem to have a ~ 1, f3 ~ 2, and £om e c 2 ~ 
a few hundred keV. From the model spectra in Figure 2, it is clear that the emission extends over 
a wide range of photon energies, from the optical to the TeV regime. We expect more optical 
emission when the dissipation takes place at larger radii. It is very likely that internal shocks 
occur over a wide range of radii, thereby extending the spectrum to longer wavelengths, with 
simultaneous optical and x-ray emission during a GRB. Recent reports of nearly simultaneous 
detection of x-ray emission from GB960720 (Piro et al. 1997), GB970815 (Smith et al. 1997), and 
GB970828 (Remillard et al. 1997) are in qualitative agreement with our expectation for internal 
shocks. 



4. Conclusions 

We have shown that the emission spectra from internal shocks are affected by synchrotron 
emission, multiple Compton scatterings, and pair creation (Fig. 1). Our model spectra could 
mimic the observed BATSE spectra (B93) for reasonable choices of the shell Lorentz factor 
7 ~ 10 2 -10 3 and shock radius r ~ 10 13 -10 14 cm (Fig. 2). 

The depletion of high energy photons due to pair creation is a sensitive probe of both T 
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and r. Detection of the high energy trough present in Figures (lb) and (2) can therefore be used 
to constrain these parameters. The potential degeneracy between the values of T and r can be 
removed on the basis of variability data, as the characteristic variability time scale depends on a 
different combination of these parameters, ~ r/2T 2 . Despite their high-energy cut-off, our model 
spectra extend all the way up to photon energies as high as 10 GeV-TeV. 

We acknowledge discussions with M. Kamionkowski, T. Piran, M. A. Ruderman, and E. 
Woods. One of us (RP) thanks R. Sari for many critical comments and stimulating discussions 
during the VIII Marcel Grossmann meeting. This research was supported in part by NASA grants 
NAG5-618 and -2859 (for RP) and NASA ATP grant NAG 5-3085 and the Harvard Milton fund 
(for AL). 
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e=hv/m e c (comoving) 



Fig. 1. — Time evolution of the radiation spectrum for r = 3.3 x 10 cm; the values of all other 
parameters are specified in the text. Panel (a) shows the spectra at different times in the early 
phase of electron cooling where only synchrotron and IC cooling are effective. The normalized 
time and radiation density A = u^{t)/uo are indicated for each curve. Because Cs/Ce = 0.33, the 
IC component begins to grow only when A exceeds ~ 0.3, as expected from equation ([l]). Panel 
(b) shows the evolution much later due to pair creation. The final density of (cold) electrons and 
positrons is n ~ 40n in this example. However, they carry only a negligible fraction (~ few percent) 
of the total energy. 
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e=hv/m e c (observed) 



Fig. 2. — Dependence of the model spectra on the dissipation radius [panel (a)] and bulk Lorentz 
factor [panel (b)] of the emitting shell. On panel (a) we show also the Band spectrum [cf. Eq. (|9|)] 
for a = 1.1, (3 = 2.15, and eom e c 2 = 200 keV, which is in a good qualitative agreement with the 
predicted spectrum for r = 3.3 x 10 13 cm. The vertical lines correspond to 10 keV and 10 MeV, 
which roughly bracket the BATSE energy range. On panel (b) we show the T-dependence of the 
spectra for a fixed r = 2 x 10 14 cm. The inferred emission radii and Lorentz factors are consistent 
with the empirical constraints set by the BATSE variability data (e.g., Woods & Loeb 1995). 



